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(N Abstract 

Q A social (sexual) net-work is modeled by an extension of the configuration model to 

Q the situation where edges have weights, e.g. reflecting the number of sex-contacts 

between the individuals. An epidemic model is defined on the network such that 
individuals are heterogeneous in terms of how susceptible and infectious they are. 
^ The basic reproduction number Rq is derived and studied for various examples, 

but also the size and probability of a major outbreak. The qualitative conclusion 
is that Rq gets larger as the community becomes more heterogeneous but that 
(~| different heterogeneities (degree distribution, weight, susceptibility and infectivity) 

can sometimes have the cumulative effect of homogenizing the community, thus 
^ making Rq smaller. The effect on the probability and final size of an outbreak is 

I— ' more complicated. 

> 

^ Keywords: basic reproduction number, heterogeneity, random network, stochastic epi- 

demic model. 

^ 1 Introduction 



X 



Epidemic models have a long history in mathematical modelling (see e.g. Diekmann and 
Heesterbeek [6] for an overview). Early models assumed a homogeneous community but 
later this was relaxed be allowing individuals to vary, for example by dividing individuals 
^ into different groups, thus defining so-called multitype epidemics. More recent models 

admitting local structures in the community have been included into epidemic models, 
household models and network models being the two main examples (e.g. Ball et al. [2] 
and Andersson llj). Admitting local structure have the effect that stochastic models are 
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favourable in that when only few individuals affect the risk of becoming infected the out- 
come should be random. 

The current paper aims at combining the two types of heterogeneities mentioned above: 
individual heterogeneities with network models. More precisely, we want to analyse how 
individual variation in susceptibility and infectivity affect the epidemic spread in a popula- 
tion composed into a social network with weighted edges. We have a sexually transmitted 
infection (STI) as a motivation for the paper. This heterogeneity between individuals 
might then correspond to varying sexual risk behaviour, e.g. not using a condom, or 
physiological features such as more or less susceptible Mucous membranes. The weights 
on the edges can for example correspond to the number of sexual contacts between the 
two individuals in question. To this end we extend the work in [5], which studies a 
weighted configuration model, which describes a stochastic epidemic model on a network, 
by including heterogeneity in susceptibility and infectivity, i.e. the ability to receive, and 
transmit, the infection. For a fixed individual i, we let the susceptibility and infectivity 
be an outcome of a pair of random variables {Xi,Yi) where these two random variables 
may depend on each other, but are assumed independent between individuals. 

We derive the basic reproduction Rq and investigate how Rq is influenced by the coef- 
ficients of variation CVx = o'x/fJ'X and CVy = ay/fiy, for different correlations, px,Y, 
between X and Y. In particular we compare the result to the Ro obtained for fixed sus- 
ceptibility and infectivity for all individuals {CVx = CVy = 0), corresponding the model 
analysed in [5]. We also investigate how the probability and size of a large outbreak 
depends on CVx and CVy- 

For an introduction to network models and their applications (including epidemics) we 
refer to Newman [11], whereas van der Hofstad [8j gives a more technical and exhaus- 
tive treaties of the random graphs and their properties, e.g. thorough treatment of the 
configuration model, and the relation between branching processes and random graphs. 
Epidemics on networks allowing individual heterogeneities have been studied earlier, but 
has perhaps not yet received enough attention. Miller [10] studies an epidemic on an 
unweighted graph where susceptibility and infectivity varies among the individuals, and 
derives bounds on the probability and size of a large outbreak. A similar problem is 
studied by Trapman and Meester in [9] who define a percolation model and derive bounds 
for e.g. expected final size and outbreak probability as function of the infectivity and 
susceptibihty. 

The rest of the paper is outlined as follows. In Section [2] we specify the graph model, the 
epidemic model and give a general expression for the basic reproduction number Rq. In 
Section |3] we investigate the effect of individual heterogeneity on Rq for different examples 
of the model. In Section |4] we derive the outbreak probability vr in a simple case and in 
Section |5] we conclude with a short discussion. 
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2 Model 



2.1 A model for weighted network 

The network model we study was originally defined by Britton et al. [5], which is an 
extension of the configuration model (see e.g. P]) but allowing for edges to have different 
weights. It is defined as follows. We have a population of size n. Individual i has a ran- 
dom number Di of half-edges where {-Di} are iid with probability function {pd(^), k>{f\ 
{Di is called the degree of individual i). Furthermore, each half-edge of individual i has a 
random weight Wij G N, j = 1, . . . ,1^2, possibly dependent of but being independent 
of each other (in the STI application Wij can for example reflect the number of sexual 
contacts individual i has with its partner j). The probability that a random half-edge of 
an individual with d half-edges has weight w, is denoted q{w\d). The network is formed 
by randomly connecting half-edges with the same weight. If the number of half-edges of 
a specific weight is odd, the last half-edge is ignored. We assume that the variance of D 
is finite. Thus the number of self loops and multiple edges are negligible if n is large (see 
Section 7 in [8]). 



2.2 An epidemic model on the weighted network 

To each individual i we assign a random vector (Xj, Yi) with probability function px,Y 
being independent of [Di.Wn, . . . ,WiD,)- The first component Xi denotes the suscep- 
tibility and Yi the infectivity (in case of getting infected) of individual i. We assume 
that < Xj < 1 and < yj < 1 and allow the variables Xi and Yi to be dependent. 
An individual i is hence described by its individual susceptibility and infectivity (Xj,yj), 
together with the independent random vector (Z)j, Wn, . . . , WjdJ indicating how many 
neighbours i has and the weights on the edges connecting to them. 



Consider two individuals directly connected to each other by an edge of weight w. Sup- 
pose the first individual has susceptibility/infectivity (xi,?/i) and that he/she is infected 
by someone else, and suppose the second individual is still susceptible having suscepti- 
bility/infectivity {x2,y2)- The epidemic model is then defined by saying that the first 
individual infects the second with probability t{w, yi, X2) defined by 

t{w,yi,X2) = l-{l-yiX2r. (2.1) 



The intuition behind Equation (2.1) is the following. In one contact, the first individual 
infects the second with probability X2yi (the more infectious the first is and the more sus- 
ceptible the second is the higher risk of disease transmission). Hence the first individual 
does not infect the second individual in w contacts with probability (1 — X2yi)^ ■ Given 
all susceptibilities, infectivities, degrees and weights, transmission events are defined to 
be mutually independent. Initially one randomly selected individual is externally infected 
(the index case) and the rest of the community is susceptible to the disease. The index 
case infects a random number of his/her neighbours following the transmission probability 
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defined in (2.1) and then becomes immune. These newly infected may in turn infect some 
of their not yet infected members and then become immune. The epidemic continues until 
there are no new infections. Then the epidemic stops. Those who were infected make up 
the final outcome of the epidemic and the number of infected is called the final size of 
the epidemic. 

We now assume that the size n of the community is large. The initial stage (before a 
non-negligible fraction have been infected) can then be approximated by a multiptype 
branching process (e.g. [7]) as is nearly always the case with epidemics in large popu- 
lations (cf. [5] for the case with unweighted edges). In the current model the type of 
the individual in the branching process approximation is characterised by the degree of 
the individual together with the susceptibility and infectivity. In order to compute the 
basic reproduction number Rq, loosely defined as the expected number of new infections 
caused by a random infected during the early stages of the epidemic, we first compute the 
expected number of individuals of a given type that an infected of a given type infects 
during the early stages of the epidemic. In branching process terminology this is the mean 
offspring matrix. 

Let nT'{di,xi,yi),{d2,x2,y2) ^e the mean number of individuals of type {d2,X2,y2) that gets 
infected by one (rfi, Xi, ?/i)-individual during the early stages of the epidemic, and let 
,yi),{d2,x2,y2)} be the mean offspring matrix. We now derive an expression 
for m(^di,xi,yi),{d2,x2,y2) ^ similar way as in [5]. Let (c?i,xi,?/i) be a fixed infected in- 
dividual in the early stage of the epidemic outbreak, and let P(d]_,xi,yi){.d2, X2,y2) denote 
the probability that it infects a (^2, 2:2, l/2)-individual along one of its di — 1 edges to 
susceptibles (the individual was infected through one of its edges). Furthermore, let 
Z{di,xi,yi){d2,X2,y2) denote the (random) number of individuals of type (^2, 2:2, 1/2) it 
infects; since type and infections along different edges are independent it follows that 
Z(di,xi,yi){d2,X2,y2) is binomially distributed with parameters di — 1 (in the early stages 
all but its "infector" are susceptible) and P{di,xi,yi){d2,X2,y2)- From this it follows that 
^{di,x^,yi),id2,x2,y2) = {d-l)p(^duxuyi){d2,X2,y2)- It remains to computc ((^2, X2, ^2) • 

We compute P(di,xi,yi){d2,X2,y2) by summing over all possible weights along the edge 
in question, since the transmission probability but also the type of the connected node 
depends on the weight. For a given weight w the probability Pw{d2,X2,y2) that the 
edge connects to a {d2, X2, ?/2)-individual is proportional to d2PD{d2) (since the fraction of 
edges connecting to (i2-individuals is d2PD{d2)/ J^d^P^id)), and proportional to q{w\d2), 
the latter being the probability that an individual of degree d2 has weight w along a 
given edge. The probability Pw{d2,X2,y2) is also proportional to px,Y{x2,y2) since the 
susceptibility and infectivity are independent of the network structure. As a consequence 
we have 

~ f, ^ q{w\d2)d2PD{d2)pX,Y{x2,y2) „^ 

Pw{d2,X2,y2) = — . I ,^ , — 7-r. • 2.2 

2^d(lMd)dpD{d) 

The quantity Pw{d2,X2,y2) is the probability that an edge with weight w connects to a 
((^2, 3^2, Z/2)-individual. We want to compute P(di,xi,j/i)('^2, a^2, 1/2), the probability that a 
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(cii, xi, ?/i)-individual infects a ((i2, 3^2, l/2)-individual along a given edge. This probability 



is obtained by summing over all possible weights, and using Pw{d2, X2, 1/2) defined in (2.2), 
but multiplied by the probability that the first individual has this weight w and multiplied 
by the transmission probability t{w,X2,yi)- We hence get 

P(di,xi,yi)(c?2,a;2,l/2) = '^qiw\di)t{w,X2,yi)PUd2,X2,y2)- (2.3) 

w 

We are now ready to compute the mean offspring matrix M with elements m(^di,xi,yi),{d2,x2,y2) 
denoting the expected number of (1^2, 2:2, |/2)-individuals that one infected {di,Xi,yi)- 
individual infects during the early stages of an outbreak. The corresponding random 
number is Z(^di,xi,yi){d2,X2,y2), and from before we know that this quantity is binomi- 



ally distributed with parameters di — 1 and P{di,xi,yi){d2,X2,y2) (defined in 2.3). It hence 
follows that 

'midi,xi,yi),{d2,X2,y2) = E[Z(^duxi,yi){d2,X2,y2)] 

= {di - l)p(di,xi,|/i)(rf2,X2,?/2) 

= (di - l)y^J{w, X2,yi)q{w\d)p^n{d2, X2,y2)- (2.4) 

The basic reproduction number, Rq, is the largest eigenvalue of the matrix M. This 
quantity plays a key role in epidemics and branching processes in that it is a threshold 
parameter. More precisely, a large epidemic outbreak is possible if and only if i?o > 1 
(see e.g. |7]). Unfortunately it is not possible to state any general features of Rq for the 
general case, except that it is the largest eigenvalue of M. In the next section we study 
some specific examples. 



3 Rq for various weight, degree, susceptibility and 
infectivity distributions 

In order to gain insight in how Rq depends on different heterogeneities: the degree distri- 
bution, the weights, and in particular the variable susceptibility and infectivity, we now 
study a few examples analytically and/or numerically. 



3.1 Unweighted network 

We begin by analysing the case with fixed weight W = 1, i.e. an unweighted network (this 
model is a special case of the model analysed in [9J). Thus we get the original configuration 
model as our social network, and the epidemic model where there is heterogeneity in terms 
of susceptibility and infectivity. In this case Equation (2.2) simplifies to 

~ . , N d2p{d2)px,Yi^2,y2) 

Pw{d2,X2,y2> = , 

where iin = '^^dpo^d) is the mean degree. 
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It follows that Equation (2.4) simplifies to 



(ai - 1)1/1 ■ (3.1) 



From (3.1) we see that the elements of M can be written as a product of two factors, 
one depending on {di,xi,yi) and the other depending on {d2,X2,y2)- The basic repro- 
duction number Rq, i.e. the largest eigenvalue of M, is derived as follows. Let {Aj} 
be the eigenvalues of M. Then trace{M) = ^2- Aj. We see that there exists vectors a 
and b (functions of {di,xi,yi), (^2, 2:2, 2/2)) such that M = a*h^ and Ma = trace{M)a. 
Thus a is a eigenvector with eigenvalue trace{M). It follows that all other eigenvalues 
must be zero and trace{M) hence equals the largest one (see e.g. j3]). To conclude, we have 



-Ro - ^ m(d,x,y),id,x,y) (3.2) 
d,x,y 

dpD{d)pxy{x,y) 



^{d-l)xy- 



d,x,y 



E{XY){E{D{D - 1)) 

2 

/ix/Xy (1 + CVxCVypx,Y) Ud + ^^^^^) . (3.3) 

\ Un / 



We have assumed that X and Y are both discrete. If they instead were continuous then 
Rq would be the largest eigenvalue of a related functional, but since we can approximate a 
continuous distribution arbitrarily well by a discrete distribution, we would get the same 



expression, so (3.3) applies whenever the network is unweighted. 



From (3.3) we see that, for an unweighted network, Rq increases with the coefficient of 
variation of the susceptibility and infectivity for the more likely scenario that infectivity 
and susceptibility are positively correlated. If on the other hand px,Y < then _Ro is 
decreasing in the coefficients of variation. 



In Figure [T] this situation is illustrated for the case D = 5, px,Y = 0.7 and px = f^y = 0.2, 
and assuming the same coefficient of variation in infectivity and susceptibility, i.e. that 
CVx = CVy. We see that Rq is increasing with CVx = CVy, which is a measure of the 
heterogeneity in the population. Larger CVx = CVy, together with positive correlation, 
results in more individuals with high ability to transmit and receive disease, and more 
individuals with low ability to transmit and receive disease. In Figure [T] it is seen that Rq 
increases with CVx = CVy in this case. Note that the distribution of {X, Y) need not be 
fully specified, only the mean and coefficient of variation of X and Y, together with the 
mutual correlation is needed. 
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Figure 1: i?o as function of CVx = CVy in the case with fixed weight W = 1 and degree 
D = 5, and where = f^y = 0-2 and X and Y have correlation 0.7. 



For the special case that X = Y, i.e. the infectivity and susceptibility are identical (fully 
correlated, px,Y = 1) but different between individuals, we get 



(3.4) 



From (3.4) we see that, in the situation where susceptibility and infectivity are identical. 



Rq increases in the randomness of the susceptibility/infectivity. 



3.2 (i-networks with random weights but fixed infectivity and 
susceptibiUty 

Suppose now that D = d and X and Y are deterministic (these are both set to ^^p, so 
that the probability of infection in one contact becomes p). The weight W is random 
with distribution q{w\d) = q{w), iid between all different edges. This is a special case of 
the model introduced in |5j. Since there is only one type of individual in this case, we get 
that 

Ro = {d-l)^t{w,,/p,,/p)q{w) 

w 

= id-l)J2a-il-pr)qiw) 

w 

= (d-i)(i-G(l-p)), 



(3.5) 



where G is the probability generating function of the weight W. Since t{w, ^/p, ^/p) = 
1 — (1 — p)™ is a concave function, we get the from Jensen's inequality that Rq = 
E{7r{W, ^Jp, ^/p)) < n{E{W), ^/p, ^/p). Thus, Rq is larger when all edges have the same 
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weight compared to the case with random weights. 



In order to obtain an exphcit expression we consider the specific example where W follows 
a negative binomial distribution W ~ NB(r,(f)), i.e. where P{W = w) = {^ZD'P^i^ ~ 
0-ju>-r^ for w = r,r + 1, . . . . The mean, variance and coefficient of variation are given by 



E{W) = fiw 



V{W) 



r(l - 



CVi 



w 



(3.6) 



1 

The probability generating function equals G{s) = (s0/(l — s(l — 0)))^. We want to 
study how Rq is affected by the randomness in the weight distribution. To this end we 
fix the mean fiw and hence set = r/fiw Having fixed fiw we see from (3.6) that 
the coefficient of variation is decreasing in the remaining parameter r. Inserting the 



probability generating function for the negative binomial distribution into (3.5) gives us 

a-p)- 



Rn 



(d-1) 1 



r 



(3.7) 



Having the transmission probability p and the expected weight fiw fixed we study how 
Rq depends on the randomness in W by how Rq depends on r. Since Rq is an increasing 
function of r (see section [6]), and CVw is a decreasing function of r, we conclude that Rq 
is a decreasing function of CVw (and in V{W)). This situation is illustrated in Figure [2] 
with = 10, D = 5 and p = 0.5, where this decay is confirmed. 




Figure 2: Rq as function of CVw where D = 5, fiw 
probability of transmission in one contact is p = 0.5. 



10 and X = r = ^05, i.e. the 



3.3 Further examples 

In the subsections above we have given two special cases of the general model where Rq 
can be computed explicitly. In general, i.e. where weights, degrees, susceptibility and 



8 



infectivity are all random, this is not the case. Below we illustrate a few "toy exam- 
ples" having all four features random, with the aim of illustrating different qualitative 
aspects of how Rq depends on the randomness of the susceptibility and infectivity as- 
suming these are positively correlated. In the examples we assume that the community 
proportion with high susceptibility is equal to the community proportion with low suscep- 
tibility and the same is assumed regarding infectivity. Varying the coefficient of variation 
in susceptibility/infectivity is achieved by varying the difference between high and low 
susceptibility /infectivity. Furthermore, the correlation between susceptibility and infec- 
tivity is managed by altering the proportion with low susceptibility and high infectivity 
(or equivalently the proportion with high susceptibility and low infectivity). In Example 
1-3 the distribution of susceptibility and infectivity {X,Y), and the degree distribution, 
is the same, but with different weight distribution. In Example 1 and 3, Rq is illustrated 
as a function of CVx = CVy, whereas in Example 2 is illustrated as a function of CVx 
with CVy fixed. 

Example 1. Our first example is where there is negative correlation between degree and 
weight (a likely scenario if degree refers to number of sex-partners and weight refers to 
number of sex-contacts per partner), and where susceptibility and infectivity are positively 
correlated. We assume that susceptibility and infectivity are random, each with two 
possible values, both having mean 0.5, such that the correlation between susceptibility and 
infectivity equals 0.8. Furthermore, the degree distribution follows a Poisson distribution 
with mean 4 truncated at 15. The weight is random with two possible values, 1 and 
10, with q{w = l\d) = 1 — d^^ = 1 — q{w = 10\d), i.e. larger degree results in higher 
probability for the smaller weight. In Figure |3j this situation is illustrated. We see that 
Rq is increasing with the coefficient of variation of infectivity and susceptibility. So, if the 
degree and weight distributions are negatively correlated, which in a sense homogenizes 
the community, then introducing heterogeneity in susceptibility and infectivity makes Rq 
increase. 

Example 2. Our second example is where degrees and weights are independent (and 
keeping positive correlation between susceptibility and infectivity). Just like in Example 
1 we let the degree follow a Poisson distribution with mean 4 truncated at 15. Also as in 
Example 1, susceptibility and infectivity each has two possible values, both having mean 
0.5, and correlation 0.8. However, in order to achieve a non monotone function of Rq, 
we let the coefficient of variation of the infectivity be fixed, CVy = 0.3. The weight is 
random with two possible values, 1 and 10, with g(l) = g(10) = 0.5, i.e. independent of 
degree. In Figure |4]i?o is plotted as a function of CVx for this situation. We see that Rq is 
not a monotone function of the coefficient of variation of the susceptibility, but increases 
initially and at one point start to decrease. This is in contrast to the situation treated in 
Example 1 where Rq increased monotonically with CVx- 

Example 3. In our last example of this section we look at the case where the degree and 
weight are positively correlated (keeping positive correlation between susceptibility and 
infectivity). Susceptibility and infectivity, as well as degree follow the same distribution 
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1.35 
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1.15 




1.1 ^ ' ' ' ' ' ' 

0.1 0.2 0.3 0.4 0.5 0.6 

Figure 3: Rq as function of CVx = CVy where the weights are random number with a 
two-point distribution. The degree D and weight W are negatively correlated, and X and 
Y are positively correlated (see Example 1 for full specification). 




2.305' ^ ^ ^ ^ ^ ^ ^ ' 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Figure 4: i?o as function of CVx-, CVy = 0.3, where the weights are random number with 
two point distribution. The degree D and weight W are uncorrelated, and X and Y are 
positively correlated (see Example 2 for full specification). 

as in Example 1. The weight is random with two possible values, 1 and 10, with q{w = 
l\d) = = 1 — q{w = 10\d), i.e. now a larger degree results in higher probability of 
the large weight. This situation is illustrated in Figure [5j We see that Rq now decreases 
with the coefficient of variation of infectivity and susceptibility. So, if individuals are 
already heterogeneous in terms of spreading the disease (degree and weight positively 
correlated), then introducing (independent) heterogeneity in susceptibility and infectivity 
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actually reduces Rq. 




0.1 0.2 0.3 0.4 0.5 0.6 



Figure 5: Rq as function of CVx = CVy where the weights are random number with two 
point distribution. The degree D and weight W are positively correlated, also X and Y 
are positively correlated (see Example 3 for full specification). 



4 Outbreak probability and size of outbreak 

The basic reproduction number Rq is not the only informative quantity of an epidemic 
outbreak. The final proportion infected, r, and the outbreak probability, vr, are also often 
of interest. For many epidemic models, including the current model, the initial phase of 
an outbreak may be approximated by a branching process. Using this approximation, the 
probability of a major outbreak coincides with the probability of the branching process 
growing beyond all limits, and this latter probability is well-known how to derive, hence 
giving a recipe how to compute vr. If, in the epidemic model, the probability for i to infect 
j is identical to the probability for j to infect i, it suffices to have one (undirected) edge 
present or not between each pair of neighbours, with the interpretation that the edge is 
present implies that if either gets infected so will the other. The resulting network of 
undirected edges then specifies who will get infected; those connected to the index case 
make out the final outbreak. Because in a large community there will be exactly one 
giant connected component (if Rq > 1, and no giant if -Rq < 1) the probability tt that 
randomly selected index case starts a major outbreak is identical to the probability that 
it belongs to the giant connected component, but the latter is simply the relative size of 
the giant connected component r. To conclude we have the well-known property that, if 
the transmission network can be modeled as an undirected network, then vr = r, i.e. the 
branching process approximation for the outbreak probability also gives the relative size 
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r of a major outbreak (e.g. Britton, 2010). 



Example 4. Consider an unweighted network [W = 1) with all individuals having the 
same degree d {D = d), and assume that individuals are of two types, with respect 
to susceptibility and infectivity, at equal proportions. Assume that both types have 
susceptibility equal to infectivity, being xi = yi = ni = n — 5 and a;2 = ^2 = At2 = + 5, 
for type 1 and type 2 respectively, with < 5 < max(/i, 1 — /i). This implies that 
t{w,yi,X2) = t{w, Xi,X2) = t{w, X2,Xi) = t{w,y2,Xi). Thus the transmission network is 
undirected in this case. Note that individuals have expected susceptibility and infectivity 
equal to fi. The offspring matrix {rriij} is defined by m^j = {d—l)fiifij/2 since an infective 
i-individual on average has {d — l)/2 susceptible type j neighbours, each of which it infects 
with probabihty fiifij {i = 1,2, j = 1,2). Simple algebra shows that the coefficient of 
variation of susceptibility/infectivity is CV = 6/ fi and Rq = /i^(l + CV^){d — 1), which 



agrees with Equation (3.4). Since the probability that the initial infective is of either type 
with equal probability 0.5, the probability of a large outbreak is vr = (tti +7r2)/2, where tTj 
is the probability of a large outbreak starting with an infective of type i. The probability 
that an i-individual infects k type 1-individuals and / type 2-individuals is 

P^{d, kJ)= (^f^ ^ (/i,/ii/2)'=(/x,^2/2)'(l - /i./ii/2 - fx,fX2/2)''-''-K 

From the theory of branching processes, we know that the probabilities tti and tt2 are 
given by 



l-TTi = 5^(l-7r;)'=(l-7r^)'Pi(^i,fc,0, (4.1) 

k+l<d 

1-7^2 = J](l-7r;)'=(l-7r^)'P2(^^,fc,0, (4.2) 



k+Kd 



where 7t[ and vTg are the solutions to the following system of equations. 



J2 {l-^[)\l-A)'Pi{d-W) (4.3) 

k+l<d-l 

1-^2 = E {^-A)\^-AfP2{d-l,k,l). (4.4) 



k+l<d-l 



The intuition behind Equations (4.3 4.4) is the following. To escape a large outbreak 
starting with one infected individual, none of the offspring may start a large outbreak. 
Thus, if an infected individual infects k type 1 individuals, and / type 2 individuals, of the 
d — 1 susceptible neighbours, no outbreak will occur with probability (1 — vr'J''(l — t^^Y . 



Taking the expected value w.r.t. the number infected of each type give Equations (4.3- 



4.4[ ). Since the index case has d susceptible neighbours (as opposed to — 1 in later 



generations) Equations (4.1 4.2) follows. 



This situation is illustrated in Figure[6]with d = 5 and /i = 0.48. We see that although Rq 
is monotonically increasing with the coefficient of variation of susceptibility/infectivity. 
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the outbreak probability vr is not monotone. Note that tt = when i?o < 1 as to be ex- 
pected. In the current example we had the transmission probability from one individual 
to another being the same in both direction. As a consequence, our transmission network 
may be modelled by an undirected network which, as argued for earlier, implies that the 
outbreak probability vr is identical to the relative size r of a major outbreak. This hence 
implies that the final size r is not monotone in the coefficient of variation as opposed to 
Rq which is monotonically increasing. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

cv 

Figure 6: The outbreak probability vr, and the basic reproduction number Rq, as function 
of the coefficient of variation of the infectiousness/ susceptibility CVx = CVy (see Example 
4 for full specification). Since the underlying transmission network is undirected vr = r, 
so the lower curve is hance also a curve of the relative size of a major outbreak if it takes 
place. 



5 Discussion 

In the current paper we have analysed an epidemic model taking place on an edge-weighted 
network, where the weights on the edges affect the transmission probability. For this model 
we have focused attention on how individual variation: in degree, weights and in particular 
susceptibility and infectivity, affect properties of the epidemic. Having no variation in the 
other factors it was seen that variation in degree increases the basic reproduction number 
Rq, and variation in susceptibility and infectivity does the same. For the situation with 
fixed degree, susceptibility and infectivity we get the opposite effect when introducing 
variation in weights. The reason for this somewhat surprising opposite effect for weight 
variation is however explained from the parametrization in the model where weights enter 
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as exponents rather than hnearly. 



An immediate hypothesis would hence be that introducing randomness in one factor when 
other factors aheady have variation, would show a similar pattern, i.e. that this would 
increase Rq- This does however not hold in general. For example it was seen in Figure [5] 
that Rq decreased with CVx = CVy when the degree and weights were positively corre- 
lated, and in Equation (3.3) when susceptibility and infectivity were negatively correlated. 
The general qualitative conclusion seems to be that Rq increases when introducing hetero- 
geneity in one aspect unless the remaining factors already induce "severe" heterogeneity. 
In a sense, if individuals are already very heterogeneous in certain aspects, introducing 
additional (independent) heterogeneity may actually have the effect that individuals be- 
come less heterogeneous (a phenomenon related to what is known as regression towards 
the mean). 



It was also observed that the probability of a major outbreak vr, and the size of an out- 
break r when the transmission probabilities were reciprocal, had a more complicated 
dependence on the randomness in degree, weight, susceptibility and infectivity. That is, 
even when Rq was monotone in CVx, the same is not necessarily true for the outbreak 
probability vr and the outbreak size r (as shown in Example 4). 

The current work may be extended in several ways in order to increase realism. For 
example, the model does not evolve in real time so the concept of length of infectious 
period and its randomness is missing. Also, if considering an epidemic over a longer time 
period, then having a time dynamic network would be of interest to study. Another feature 
of interest would be to study the effect of vaccination, or other preventive measures, on 
the reproduction number. Further work might also analyse the current model in more 
detail and/or construct other examples, to further investigate the qualitative properties 
of the model. 
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6 Appendix 



Let c = yz^, then it follows from Equation (3.7) that 



i?o = (rf-l)(l- (- 



il-p) — 



I — p)— + p 

[d — 1)(1 — exp{^—rlog{\ + c/r)}) 

d{rlog{\ + c/r) 



(i?o)'r. = -exp{-rlog{\\ clr)}{^ 

= exp{—rlog{l + c/r)}(log{l + c/r) 



c/r 



1 + c/r 



for all r > 0. Define z = c/r and g{z) = log{l + z) — where 2; > 0. It follows that 
g'{z) = 1/(1 + z)- 1/(1 + zf. Since ^(0) = g'{0) = and g'{z) > for z > 0, it follows 
that g{z) > when 2; > 0. We conclude that (-Ro)r > when r > 0. Thus Rq is an 
increasing function of r. 
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